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Abstract 

This is the hrst of a series of papers in which a new formulation of quantum 
theory is developed for totally constrained systems, that is, canonical systems in 
which the hamiltonian is written as a linear combination of constraints ha with 
arbitrary coefficients. The main purpose of the present paper is to make clear that 
classical dynamics of a totally constrained system is nothing but the foliation of 
the constraint submanifold in phase space by the involutive system of inhnitesimal 
canonical transformations Ya generated by the constraint functions. From this 
point of view it is shown that statistical dynamics for an ensemble of a totally 
constrained system can be formulated in terms of a relative distribution function 
without gauge hxing or reduction. There the key role is played by the fact that the 
canonical measure in phase space and the vector helds Ya induce natural conserva¬ 
tive measures on acausal submanifolds, which are submanifolds transversal to the 
dynamical foliation. Further it is shown that the structure coefficients dehned 
by {ha,h^} = J2-yCaph^ should weakly commute with ha, J2'Y{h-y,Ca^} ~ 0, in or¬ 
der that the description in terms of the relative distribution function is consistent. 
The overall picture on the classical dynamics given in this paper provides the basic 
motivation for the quantum formulation developed in the subsequent papers. 


1 Introduction 



In the canonical approach to quantnm gravity the Dirac quantization prescription oJ 
is widely adopted, in which the classical constraints ha ^ ^ are formulated as con¬ 
straints on physical states of the form 


( 1 . 1 ) 


ha^ = 0 . 


As is well-know, however, this leads to the frozen formalism in which the dynam¬ 
ical evolution equation is lost. This problem is closely related with the fact that 


operators corresponding to the physical quantities which play the role of time are 
excluded from observables in this . Further the Dirac quantization 



of general relativity does not give a mathematically well-dehned formulation apart 


from the regularization and the operator ordering problems because Eq. (|l.lj) does 


not have normalizable solutions in general even if the spatial diffeomorphism free¬ 
dom is eliminated before quantization. 

As discussed in Ref. [@ , this difficulty comes from the too formal application of 
the Dirac procedure. Since the hamiltonian is written as a linear combination of 
the constraint functions in general relativity, the hamiltonian constraints carry all 
the information on dynamics. Hence if we formulate the hamiltonian constraint 
as the condition on state vectors as above, each state vector becomes a dynamical 
object. This should be compared with the ordinary quantum mechanics. There 
state vectors are used to describe dynamics of a system, but each state is not a 
dynamical object. Dynamics is described by a one-parameter family of states <I>(t) 
satisfying the Schrodinger equation. Each state vector in this family merely carries 
information by possible measurements at each instant t. 

This observation indicates that we should impose the hamiltonian constraints 
on the object which picks up all the possible state vectors allowed by dynamics. The 
most natural object of this nature will be the probability amplitude \h($) which 
assigns the probability to each state vector, taking account of the probabilistic 
nature of quantum mechanics. 

Here one may notice that a similar phenomenon occurs in the classical dynamics 
of a system with hamiltonian constraints. For simplicity let us consider a system 
with a single hamiltonian constraint h 0 on a phase space F. If we reduce this 
system to a canonical system without constraint, dynamics is described by a curve 
7 o(f) in a reduced phase space Fq which is a solution to a canonical equation of 
motion. This curve corresponds to the family of states <h(t) in quantum mechanics. 
On the other hand in the original phase space this curve corresponds to a curve 7 
contained in the constraint hypersurface which corresponds to T above. This 
analogy becomes better if we consider an ensemble of systems instead of a single 
system, for which 70 (t) is replaced by a family of distribution functions po(t) on 
To and 7 by a distribution function p on F. Clearly p does not represent a state 
but is a dynamical object which picks up possible states allowed by dynamics. 
Hence p should be constant along each curve in F corresponding to a solution 
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to the equation of motion in Fq. This implies that it is not normalizable on F. 
This may be regarded as the essential reason why the solutions to Eq . (11.11) is 
unnormalizable. It is obviously absurd to postulate dynamics so that it picks up 
p which is normalizable in F. Applying the quantum hamiltonian constraints on 
states just corresponds to such an approach. 

This observation suggests that the investigation of the structure of classical 
statistical dynamics of totally constrained systems will shed a good light on how 
to hnd a consistent quantum formulation of them. This is the motivation of the 
present paper. Since a variety of forms exist for the canonical formulation of gravity 
and since the structure of the problem is common in all the theories with general 
covariance, we consider a generic totally constrained system in most part of the 
paper. 

The organization of the paper is as follows. First in §2 we consider a simple 
totally constrained system obtained by embedding a canonical system without con¬ 
straint to a larger phase space in order to find how to interpret the unnormalizable 
distribution function on the extended phase space. Then in §3 on the basis of the 
result obtained there we describe how to formulate the statistical dynamics of a 
generic totally constrained system with a single constraint without reduction and 
not referring to special time variables. Further the general structure of reduction 
and its freedom is examined because it is relevant to the time variable problem. In 
particular by applying it to the totally constraint system describing a relativistic 
particle in curved spacetime, it is shown that the background spacetime should 
have a Killing vector in order that there is a natural reduction of this system. In 
the subsequent two sections the formulation obtained for a single constraint system 
is extended to a multiple constraint system. First in §4 an overview on the canon¬ 
ical structure of general relativity in terms of the ADM variables is given in order 
to make clear that dynamics of a totally constrained system with multiple con¬ 
straints is completely determined by the foliation of the constraint submanifold by 
the involutive system of the inhnitesimal canonical transformations generated by 
the constraint functions. Then in §5 based on this viewpoint, statistical dynamics 
for multiple constraint systems is formulated in terms of the relative distribution 
function by proving the existence of natural conservative induced measures on 
acausal submanifolds. Section 6 is devoted to discussion. 


2 Embedding of an unconstrained system into a 
constrained system 

In this section we study the dynamics of a simple constrained system obtained by 
embedding an unconstrained canonical system into a larger phase space. The main 
purpose is to hnd the way how to formulate the dynamics of a constrained system 
and its ensemble without reducing it into an unconstrained one. 
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2.1 Canonical system 

A canonical dynamical system with no constraint is specified by a triplet (r,a;, h) 
of a phase space, a symplectic form and a hamiltonianB. The phase space T is a 
2n-dimensional smooth manifold, the symplectic form a; is a closed non-degenerate 
2-form on T, 


duj = 0, (2.1) 

;= i-A*" 0^ = ^ 0, (2.2) 

n\ n\ 

and the hamiltonian h is a smooth fnnction on x T. 

Let JT(r) and r¥(r) be the sets of all smooth fnnctions and all smooth vec¬ 
tor fields on T, respectively. Then for any / G .?^(r), uj nniqnely determines a 
vector field Xf G A’(r), which is called the infinitesimal canonical transformation 
generated by /, throngh the eqnation 

df = -IxfUJ ^ Vf = u;{V,Xf) Vl/GA(r), (2.3) 

where ly is the inner prodnct operator which maps a p-form y; to a (p — l)-form 
defined by 

(IxxWi, • • •, Vi) = 14, • • •, Vi)- (2.4) 

From the identity 

X ~ Ix o d -|- d o (2-5) 

the infinitesimal transformation Xf preserves uj and fl; 

£;XfJJ — 0 , iQxf^ ~ 0 . ( 2 - 6 ) 

Conversely any vector field which satisfies this eqnation is an infinitesimal canonical 
transformation generated by some fnnction on T, at least locally. 

In terms of this infinitesimal canonical transformation the Poisson bracket of 
two fnnctions / and g are defined by 

{/,9}:=-A> = -a,(.Y,,Ag. (2.7) 

It follows from this definition that 

•V/.») = -lV/.AJ. (2.8) 

Thns the correspondence f ^ Xf gives a homomorphism from JF(r) into A’(r) as 
Lie algebras, whose kernel is constant fnnctions. 

In appropriate local coordinates (qfipj) uJ can be always written as 


UJ = dpj A dqfi 


( 2 . 9 ) 
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( 2 . 10 ) 


In this local coordinate system Xj is expressed as 

X 

dpj dq^ dq^ dpj ’ 

which leads to the familiar expression 

If 1 df dg df dg 

’ ^ dq^ dpj dpj dq^ ’ 


( 2 . 11 ) 


Finally the dynamics in the phase space is determined by the hamiltonian h 
in the following way. Let the canonical coordinate of i? in J? x F be f, and for 
each value of t let Y := be the infinitesimal canonical transformation generated 
by h regarded as a function on F. Then for a single system, its possible histories 
are given by the integration curves of the vector held dt + Y on R x T when t is 
regarded as the time variable, and the value of / G JF(F) along each curve, when 
regarded as a function of t, satishes the equation 


In particular in a local coordinate system [t, u°‘) of i2x F each integral curve follows 
the canonical equation of motion 


du°' 

dt 


= {u^,h}. 


(2.13) 


On the other hand the behavior of an ensemble of the system is described by a 
distribution function p G tF{R x F), which satishes the two conditions: i)/p pVl = 1 
for each f, and ii) pVl is preserved by the equation of motion. This second condition 
yields the equation of motion for p. 


;ea,+v(pO)=0 (a, + F)p = 0, (2.14) 


which implies that p is constant along each integration curves in x F. When 
we introduce pt dehned by ptiu) = p(t,u){u G F), the expectation value of / G 
X{R X F) at a time t is given by 



< / >t= fpM, 

(2.15) 

and its time-derivative by 



d . 

^ L ^ It (f 

(2.16) 



(2.17) 


In particular, when / is a constant of motion, i.e., -£ = 0, < f >t becomes time- 
independent. 
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2.2 Embedding into a totally constrained system 

The above canonical dynamical system can be embedded into totally constrained 
canonical systems with larger phase spaces by varions ways. Here we consider the 
simplest one. 

Let (r, d), h) be a canonical system defined by 


T ■.= E? xT 3 {q^,po,u), 

(2.18) 

Cj := dpo A dq^ -f oj, 

(2.19) 

Hq°:Po,u) := h{q°,u)+po. 

(2.20) 

Then the infinitesimal canonical transformation F generated by h 

is expressed as 

Y = Y + dgO- dgohdp,. 

(2.21) 

Hence by the projection vr defined by 


vr : f R X T , 

(2.22) 

LU LU 


{q°,Po,u) (g°,n) 



it is mapped to = Y + dt, and each integral cnrve of F to a solntion to 
the eqnation of motion in i? x T. Therefore, noting that h is conserved and vr is 
injective on each h =const snrface, one sees that the original canonical system is 
eqnivalent to the extended canonical system with a constraint h =const. 

In this embedding only the integral curves of F in the extended phase space, 
which we call the hamiltonian flow, have a physical significance, while the canonical 
time variable for the extended system, which is denoted by r, just plays the role 
of a parameter of these curves. Hence for an arbitrary function N{t) the system 
with h replaced by if = Nh also equivalent to the original system under the 
constraint h =const. In particular for the special choice of the constraint, h = 0, 
this equivalence holds for an arbitrary function N G J^{R x T) since 

= A^F + IXn ~ NY. (2.23) 

where A k, B means that A = B under the constraint. We express this last 
situation by saying that the original canonical system is embedded into the totally 
constraint system (r,d;, h). In this expression h is understood to play the double 
roles, one as the generating function of the hamiltonian flow and the other giving 
the constraint h = 0. 

2.3 Distribution function on the extended phase space 

The distribution function p for the unconstrained system, if it is regarded as a func¬ 
tion on T, is constant along the hamiltonian flow from Eq. (p.l4|) and Eq.( p.2l|) . 
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This is quite natural since each pure dynamical state of the constrained system is 
not represented by a point but by a hamiltonian flow line in the extended phase 
space. Hence, taking account of the constraint, it is natural to consider the distri¬ 
bution function p on the extended phase space defined by 

p = p5{h), (2.24) 

where 5(*) is the delta function. From this definition it follows that p is charac¬ 
terized as a distribution on F satisfying the two equations, 

Yp = 0, (2.25) 

hp = 0. (2.26) 

Here by a distribution on F we mean a functional / on the space of smooth functions 
(j) with compact supports in F, which is expressed as 


<f,4>>-= pm 


(2.27) 


when it can be identified with a function on F. 

Since the original phase space F at time t can be identified with the intersec¬ 
tion of the = t hypersurface St and the constraint hypersurface Hh in T, the 
expectation value of f G J^iRxT) at a time t is expressed in terms of p and 
Et := 5{q^ - t) as 

<f>t=<EtpJ>. (2.28) 


This fixes the interpretation and the normalization of the distribution function p 
in the extended phase space. 

From the dynamical point of view these =const surfaces have no special 
significance in the extended phase space, apart from that they are ‘Cauchy surfaces’ 
for p. In fact we can easily extend the expression |2.28| to that for the expectation 
value on an arbitrary hypersurface S which is transversal to all the hamiltonian 
flow lines. Let us call such a surface a maximal acausal hypersurface and denote 
the corresponding expectation value by < / >s. Then for two maximal acausal 
hypersurfaces Si and S 2 , < / >Si and < / >S 2 should coincide for any constant 
of motion /. If we express < / >s in terms of a measure ps as 


< / >s— 



(2.29) 


the above condition implies that = /^S 2 where 9 is the diffeomorphism from 

Si onto S 2 determined by the hamiltonian flow. Hence the requirement that 
coincides with that given by Eq. (p.28|) completely determines /is for any maximal 
acausal hypersurface. 

The explicit form of the measure is given by 




(2.30) 
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where |x|s implies that the differential form x is regarded as a positive measure 
on S. In order to see this, first note that, for any constant of motion /, //phi is a 
closed form from 

difiyCi) = £p(/hi) - lydifn) = 0. (2.31) 

Hence by applying the Stokes’ theorem to a region bounded by Si and S 2 one gets 
fdiJ.Y ,2 ■ On the other hand from the identity 

dcf) A Jyhi = —Iy{d(l) A hi) + lydcpi} = —{h, 0}O, (2.32) 

one can rewrite the above expression for < / >s as 

< / >s=< Ej:P, f >; 

/;s= |{h,0}|h(0-r), (2.33) 

if the maximal acausal surface S is specified by the condition 0 = r(=const). It 
is easy to see that the right-hand-side of this equation coincides with Eq.( |2.28|) for 
0 = and t = t. 

2.4 Probability interpretation of the relative distribution 
function 

Since p is a distribution on F, we are tempted to interpret it simply as giving a 
probability measure for measurements of physical quantities defined on F. Follow¬ 
ing this interpretation, the probability for a set of quantities fi, ■ ■ ■ E E(T) to take 
value in a given set of ranges Ai, ■ ■ ■ C it should be given by 

Pr(/i eAy-) = C<p,Ef^ (Ai) ■ ■ • >, (2.34) 

where C* is a normalization constant and Ef{A) is the characteristic function of the 
region where the value of / is contained in A. However, this naive interpretation 
does not work by itself because the integration of p over the whole extended phase 
space diverges: 

< p, 1 >= J dt < 1 >t= 4-00. (2.35) 

Hence p cannot give a hnite measure on F by simple renormalization. Nevertheless 
it can be interpreted as giving a relative probability density under some limited 
situations. 

To see this, let us dehne the conditional probability for physical quantities 
/i, • • ■ to take value in Ai, • • ■ under the condition 0 = r, by 

Pr(/i £ Ai, • • ■, 0 G A^) 

A 0 e ^ 4 ) 

T e A^ 
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Pr(/i G Ai, • • ■ 10 = r) : 


(2.36) 





Then the expectation value of any / G JT(r) determined from this probability 
coincides with the right-hand side of Eg. (|2.33|) , if and only if {h, 0} is given by 
some function of 0 on T^h- This condition requires that 0 is expressed in terms of 
a function g and some constant of motion k as 

+ (2.37) 

Let us call a function satisfying this condition a good time variable. 

When a maximal acausal hypersurface S is given, we can always find a good 
time variable which is constant on S. However, its freedom is just a rescaling of 
the variable so that it is in general impossible to find a good time variable which is 
constant on each of more than two given acausal hypersurfaces. Though this fact 
has no significance in the classical framework, it seems to have a deep implication 
in the quantum framework in connection with the unitarity problem as will be 
discussed in the next paper®. 


3 Non-trivial system with a single hamiltonian 
constraint 

In this section we show that by a simple generalization of the formulae in the 
previous section we can discuss the dynamics of a totally constrained system with 
a single hamiltonian constraint without reducing it to an unconstrained system. 

3.1 Dynamics in the extended phase space 

Let (T, 00 , h) be a totally constrained system with a single hamiltonian constraint 
h = 0, and Y G r¥(r) be the infinitesimal canonical transformation X^. Then 
from the consideration in the previous section, it is natural to interpret that each 
integration curve of Y on the constraint hypersurface Hh yields a possible time 
evolution of the system. Hence if Y has a zero- point on it represents a solution 
for which any physical quantity takes a fixed value. Since such a solution is quite 
unphysical, we assume that Y does not vanish on On the other hand if there 
is a closed orbit in the hamiltonian flow, it represents a completely periodic world 
like the anti-de Sitter spacetime. We do not consider such causality violating cases 
in this paper either. We further assume that the hamiltonian flow is not ergodic. 
This is equivalent to requiring that the foliation of by the hamiltonian flow has 
a locally trivial bundle structure whose fiber is homeomorphic to R. Thus it has a 
global section. We call an extension of such a global section off T^h as hypersurface, 
a maximal acausal hypersurface. 

When we consider an ensemble of totally constrained systems with the same 
structure, they are represented by a set of hamiltonian flow lines in T, each of 
which intersects with a maximal acausal hypersurface at a single point. In the 
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limit that the ensemble consists of very large nnmber of members, these intersection 
points determine a measure ps on each maximal acausal hypersurface S. From 
its dehnition this measure is preserved by the mapping among maximal acausal 
hypersurfaces determined by the hamiltonian flow. Let c be a positive dehnite 
constant of motion. Then, since c/yff is a closed form and yields a measure with 
the same property on each maximal acausal hypersurface as shown in the previous 
section, the Radon-Nykodim derivative of ps by c|/yff| yields a function which 
is constant along the hamiltonian flow on the constraint surface Thus if we 
regard this function as a distribution whose support is contained in E//, we are 
naturally led to the distribution function p on T which satishes 

Fp = 0, (3.1) 

hp = 0 (3.2) 


in the distribution sense. From its dehnition the expectation value of a physical 
quantity / G JF(F) on a maximal acausal hypersurface E for the ensemble is given 
by 

< / >s= [ /cp|/yfi|/ [ cp|/yfl|. (3.3) 

In practical situations each maximal acausal hypersurface is specihed by the 
condition (j) =const in terms of a physical quantity (j) on the phase space F. Here 
the constant should be a special value and the other values may not give maximal 
acausal hypersurface in general. We call such a function instant function. 

In terms of the instant function the dynamics of the totally constrained system 
is formulated in the following way. First one selects an appropriate instant func¬ 
tion 01 from measured quantities, and specihes a maximal acausal hypersurface, 
say, by 01 = 0. Then measurements of various physical quantities determines a 
measure psj, which in turn determines the value of the distribution function p 
on the maximal acausal hypersurface Ei. p is uniquely extended over the phase 
space by Eqs.( p.l|) -( p.2|) , at least around the constraint hypersurface Once the 
distribution function is determined, one can calculate the expectation value of any 
physical quantity at an instant specihed by any instant function. Here, though 
the distribution function depends on the choice of the constant of motion c in the 
above procedure, this freedom does not affect the predictions on the expectation 
values since c and p come into the theory always in the combination cp. 

The reason why we have introduced the apparently superhuous freedom of c 
in the dehnition of p is to widen the concept of good time variables introduced 
in the previous section. Let 0 be an instant function such that 0 = r(=const) 
gives a maximally acausal hypersurface for any r in some open interval of R. If 
we require that the probability measures on these hypersurfaces derived from the 
natural measure pQ as in Eq. (p.36|) coincides with the conserved measure ps, we 
obtain the condition 

c{h,0} /(0), (3.4) 
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where /(0) is an appropriate function of 0. This condition is equivalent to the 
condition that 0 is a function of 0o which is a solution to the equation 

cy0o~ const. (3.5) 

We will show later that we must allow for a nontrivial choice of c for a natural 
time variable to satisfy this condition even in simple cases. 

3.2 Reduction 

As clarihed in the previous section, the dynamics of a totally constrained system 
{r,u,h) can be described with no reference to a special time variable. Now let 
us study the relation of this description with that in terms of a time variable in a 
canonical system with no constraint which is obtained by reduction. 

In general a canonical dynamical system without constraint (ro,C(;o,ho) is a 
reduction of the totally constrained system (r,C(;, h) if there exists a diffeomorphic 
embedding 

$ ; R X Tq —^ S/f (3-6) 

which satisfies the following two conditions: 

i) ^M + yo) = kY {keJ^{T)), 

ii) Z 2 ) = ujq{Zi, Z 2 ) Zi, Z 2 G r¥(f X To), 

where Yq is the inhnitesimal canonical transformation on Tq generated by /iq. A 
convenient characterization of <h is given by the following well-known result. 

Proposition 3.1 The necessary and sufficient condition for the mapping <h to give 
a reduction is that the following equation holds: 


<h*a; = c^o — dhf) A dt. 


(3.7) 


Proof 

1) Necessity. 

From the condition ii) for the reduction mapping there exists a 1-form ^ such 
that <F*a; = Uq — f A dt. Since u and ci^o are both closed, we obtain df A dt = 0. 
From this it follows that we can choose f so that it is closed. Hence <F*C(; can be 
written as <F*a; = Uq — dp A dt in terms of some function p. Applying Io^+Yo on this 
expression leads to 


l3^+Yo^*uj = d{p - ho) - {dtP + Yop)dt. 
On the other hand from the condition ii) 

Iat+Yo^*oj = ^*{klYUj) = -^*{kdh) = 0. 
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Hence we obtain dp = dho + {dtp + YQp)dt, which leads to the equation in the 
proposition. 

2 ) Sufficiency. 

If the equation in the proposition holds, the condition ii) is obvious. Further it 
is easy to check that <h*/<i>^(a 4 +yp)Ci; = lot+Yo^*^ = 0. Hence there is a function k 
on r such that 

I^*{dt+Yo)^ ~ —kdh = IkY^- 

This is equivalent to the condition i). I 

The reduction mapping $ induces a function 0 and a vector held V on Ejj- 
dehned by 

t = $,9* = H, (3.8) 

which are related by 

V(j) = 1. (3.9) 

The function 0, when extended off the constraint surface J^h, yields a time variable 
on T, and foliates T into a family of acausal hypersurfaces = {m G r|0(M) = t}. 
On the other hand the vector held V, when extended oh generates a one- 
parameter family of transformations on T such that = ^t+r, and JF(ro) 

can be identihed with the restriction to T^h of the set of functions on T that are 
invariant by these transformations. Hence the reduced phase space Tq can be 
naturally identihed with SoUSiif. Under this identihcation the reduction mapping 
$ can be written as 


$ ; RxVq - 

■ 

(3.10) 

If 

If 


{t,x) 




Hence, when an acausal hypersurface Sq is given, the vector held V completely 
determines the reduction. For this reason we call V the reduction field. 

The reduction held yields a reference dynamics in describing the dynamics of the 
totally constrained system by time evolution, and the hamiltonian in the reduced 
system is essentially the generator of the deviation of the hamiltonian how from this 
reference dynamics. To see this, let 9t be the mapping from Sq to St determined 
by the hamiltonian how. Then pt = yields a family of transformations on 

To = So n S/t describing the deviation of the hamiltonian how and the reference 
dynamics. For a given point m G F the tangent vector Z^pu)if) of the curve ptiu) 
is given by 

— ~Yrn{u) + = (z/_t)*( —U-|- kY)Qp.Y)i (3-11) 

where /c is a function determined by the condition kY k , 1 . Hence from the condi¬ 
tion ii) for the reduction mapping we see that Z{t) is nothing but the inhnitesimal 
canonical transformation Tq generated by h^. 

An arbitrary vector cannot be a reduction held. It must be approximately a 
canonical transformation as the following proposition shows. 
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Proposition 3.2 A vector field V onT is a reduction field if and only if it satisfies 
the following condition: 

i) Vh ^ 0, 

a) there exists a function f on T such that Vf ~ 1 and £ya; A = 0. 

In particular any infinitesimal canonical transformation satisfying i) is a reduction 
field. For this case the reduced hamiltonian h^ can he chosen to he independent of 
t. 

Proof 

1) Necessity 

From the definition of the reduction held i) is obvious, and for the time variable 
(j) in the reduction the former half of condition ii) is satished by dehnition. Further 
since <F*a; = ojq — dh^ A dt from Proposition 

^*{dlvuj A df) = dIdt^*oj A dt = 0. 

This is equivalent to the latter half of ii). 

2) Sufficiency 

By taking the function in condition ii) as a time variable, let us construct the 
mapping <F : x Fq — Hh from and V exactly in the same way as described 

below Eq. (|3.9|) . Then it obviously follows that = V and <F*0 = t. Hence from 
ii) we obtain 0 = ^*{dIvoj A df) = d{Id^^*uj) A dt. This implies that there exist 
hQ,a E IF{R X Fq) such that = dh^ + adt. Since Jl- = 0, a is expressed as 

a = —dthQ. Hence <h*a; is written in terms of 2-form ujq such that Io^Uq = 0 as 


<h*a; = cuo — dho A dt. 

Since oj is closed, coo is also closed. Further from 

£:q^ujq a dt = A dt = A dcfl) = 0, 

we obtain = 0. Hence from Proposition $ yields a reduction mapping. I 

Though the reduction helds are rather restricted, we can hnd a reduction $ for 
which ho = <F*p, for an arbitrarily given function p on P. To see this, let 0 be a 
function on P such that Y f = —{h, 0 } 7 ^ 0 on J2h, and put V = Xp + aX^ where 
p and a are functions to be determined so that H is a reduction held. First from 
condition i) in Proposition a is uniquely determined as a = —{h,p}/{h, 0}. 
Since = —da A df, condition ii) simply reduces to Vfi = —X^p Ri 1. As 

X^ is transversal to 12h, we can always hnd a solution p to this equation such 
that p coincides with an arbitrarily given function on J2h. Since ho = and 

Vp ~ Yp/Y(j), this implies that ho can be an arbitrary function on i? x Pq. 

This observation shows that reduction and the corresponding reduced hamil¬ 
tonian has a physical signihcance only when the original system has some kind of 


13 


time-translation symmetry which induces a reduction held in the phase space. In 
generic cases for which no such symmetry exits it is more natural to discuss the 
dynamics in the original phase space with a constraint as described in the previous 
section. 


3.3 Example: a relativistic particle in curved background 

We illustrate the argument so far by a simple totally constrained system describing 
a relativistic particle moving on a curved background (M, g) with a mass which 
may be position-dependent. 

The action for this system is given by 

82 = -2 J dT[-U{x)g^^{x)x>"x'"]^^‘^. (3.12) 


This action is equivalent to 

Si = J dr {p^{x^ -v>^) -2 , (3.13) 

where is a Lagrange multiplier. The variation of this action with respect to 
yields 

P^. = ^ ^ . (3.14) 


By eliminating in Si with the help of this equation we get 

1 


S = J dr Pf,x^ - -N{g^’^p^p^ + U) 


(3.15) 


where N is regarded as an independent variable. Hence the original system is 
equivalent to the totally constrained canonical system 


(r,o;,h); r = T*M, 

uj = dO = dp^ A dx^, 

h = ^ [g^''{x)p^p^ + U (x)], (3.16) 

where T*M is the cotangent bundle of M and 6 is its canonical 1-form. 

This totally constrained system is not simply reducible as the hamiltonian con¬ 
straint is quadratic in the momentum unlike the system considered in §2. However, 
one can discuss its dynamics in the sense discussed in this section except for special 
cases. In fact, the generating vector Y of the hamiltonian flow of this system is 
given by 

„ (9 1 / dg^"^ dU \ d , 
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which vanishes at points where = 0 and dU/dx^ = 0. Hence if there is no 

point such that U = 0 and dU/dx^ = 0, the hamiltonian flow has global acausal 
hypersurfaces. 

As shown in the previous section, we can always find a reduction of this system 
into an unconstrained system. However, in order for the corresponding reduction 
field to be associated with some time translation of the system, the system must 
have a symmetry. To see this, let A' be a vector field on M. Then it induces a 
vector field K on the phase space T*M which is expressed in the local coordinate 
system {x^,pi,) as 

It is easy to see that this is an infinitesimal canonical transformation generated by 
the function p^K^. The condition for this field to be a reduction field is given by 
the following proposition. 


Proposition 3.3 The vector field K yields a reduction field of a totally constrained 
system {T*M,u!, h) if and only if Kh ~ 0. In particular for the constraint h given 
by Eg. (\3. 1 6[j this is eguivalent to the condition that K is a Killing vector of the 
metric g = Ug. 


Proof 

Since K is an infinitesimal canonical transformation the former half of the 
proposition is obvious from Proposition |3]^. For h given by Eq. (|3.16|) the condition 
is expressed as 


0 ^ k{U-^h) = -V^{g,^K^)~g^^~g^ 


PxPo 


where is a covariant derivative with respect to g. This equation is equivalent 
to 


= 0 , 


which implies that iP is a Killing vector of the metric g. I 


This restriction on the system is closely related with the condition for the 
system to have a good time variable which is independent of p^. To see this, let cf 
be a function on M. Then, since Ycf is now given by 


the condition for 0 to be a good time variable, Eq. (|3.5|) , is written as 

^^ 0, (3.20) 

which is equivalent to 

(3.21) 
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This last equation is written in terms of the covariant derivative with respect to 
the metric as 

+ V,(17V^0) = 0. (3.22) 

Hence the metric g^i, must have a Killing vector K such that 

K^^ = Vy. (3.23) 


In particular must be static. 

Putting these two requirements together, we hnd that there should exists a 
function 0 such that is a Killing vector of g and =const in order 

that there is a good time variable which is a function on M and whose gradient 
held generates a reduction held. It is easy to see that these conditions are satished 
if and only if g has a static Killing vector along which U is constant (cf. Kuchaf’s 
article in Ref. |^. For such cases the reduced hamiltonian is given by 

For example for a relativistic free particle in Minkowski spacetime for which 
U = =const, the translation Killing vector satishes these conditions when 
is a constant time-like vector. The time function and the reduced hamiltonian 
are given by 0 = and ho = a^p^- On the other hand the boost and rotation 
Killing vectors do not correspond to good time variables on M because they are 
not gradient vectors. 


4 General relativity as a totally constrained sys¬ 
tem 

Before extending the argument on the totally constrained system with a single 
constraint to a more generic case, we give an overview on the structure of the totally 
constrained system obtained from general relativity The main purpose is to 
make clear that the classical dynamics of general relativity as a totally constrained 
system is nothing but a foliation of the constraint submanifold such that each leaf is 
one-to-one correspondence with a 4-dimensional diffeomorphism class of solutions 
to the Einstein equations. This fact will be used to establish an interpretation of 
generic systems in the next section. 

4.1 ADM canonical formulation on the 3-metric space 

For simplicity we only consider globally hyperbolic vacuum spacetime {M,g), and 
assume that it is spatially compact. Hence M is diffeomorphic to it x S' where S 
is a compact space. Let q(t) be the induced 3-metric on {t} x S, K(t) the extrinsic 
curvature tensor of {f} x S, and n(t) = — v) the unit normal to {t} x S where 

p G r¥({t} X S). Then by regarding g(t), K{t), N{t) and vit) as quantities on S 
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with the time parameter t, K{t) is written as 

^jk 2,N ^ Djl^k (4-1) 

where Dj is the covariant derivative with respect to q. Further the Einstein-Hilbert 
action for the spacetime {M,g), 

is expressed in terms of these quantities as 

l^d^x^\N{m + g{K, K) - (IrKf). (4.3) 

By introducing the momentum variable conjugate to qjk by 

^k ^ _ gi^TriF), (4.4) 

ZiKj 

This action is put into the canonical form, 

S = J dt[< p, q > — {< hj_, N > + < hr), v >)], (4.5) 

where h±, hj) and p are linear functionals of functions, vector helds and 2nd-rank 
covariant tensor fields on S, respectively, dehned by 


< h±J >:= J^d^xf (^g(p,p) - ^(Trp)- 

< ho, X >:= d^x{-2XWkp)) =< p, &xq >, 
<p,v >:= f d^xvjkp’^. 



(4.6) 

(4.7) 

(4.8) 


Let T^{S) be the set of all smooth (p, g)-type tensor held on S', and g be a 
reference Riemannian metric on S. Then by taking the completion with respect to 
the inner product 

{a,(3) = ^ j d^x^\q{b^a,&(3), (4.9) 

we obtain the Sobolev space H\^(Tp{S)) which is a real Hilbert space. In particular, 
if we dehne the space of 3-metrics on S by 

AI(S') := |g G i7^(7^°(S'))|g is positive dehnite on S'! , (4-10) 
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A4{S) becomes an open snbset of H^iT^^S)). Hence it has a natnral Hilbert 
manifold strnctnre, and its tangent space and cotangent space are both isomorphic 
to the original Hilbert space: 


T,iM{S)) - T;{M{S)) - Hl{T,\S)). (4.11) 

Farther, since the operator L dehned by 

{a, (3) =< a, LP >, (4-12) 

where = q^^q’^'^Pim-, is given by L = and elliptical, it dehnes an 

injection from T*(M(S)) into the space of linear fnnctionals on Under 

the identihcation by this mapping the momentnm p can be regarded as an element 
ofT;(M(^)). 

From this observation the action Eg. ([4.51) dehnes a canonical system (F,a;,iJ) 
with an inhnite nnmber of constraints: 

< ho, X>=0 VX e (4.13) 

<h^J>=0 V/eX(^), (4.14) 

where the total phase space F is given by the cotangent bnndle T*PA{S), the 
symplectic form formally by uj =< 6p A 6q >, and the hamiltonian H by 

H =< h±, N > + < V > . (4-15) 

Thongh we can give an exact expression for uj by introdncing the basis of H^(TP(S)), 
we will not do it becanse the argnment in this section is formal. Since the hamilto¬ 
nian is written as a linear combination of the constraint fnnctionals, this canonical 
system is a totally constrained system. 

4.2 Dynamical foliation of the phase space and Diffo(M)- 


classes 

Let E D and be the snbmanifolds of F dehned by 

Ed := {n e F| < X > (u) = OVX e 4’(^)} , (4.16) 

Eh := {n e F| < / > (n) = OV/ G X(^)} . (4.17) 

Then from the Poisson bracket strnctnre among the constraints, 

{< ho, Xi >,< hn, X 2 >} =< ho, [Xi, X 2 ] >, (4-18) 

{< hD, X >,< hj_, / >} =< h_i_, £rxf >, (4.19) 

{< hD, /i >, < hD, f2 >} =< hD, fiDf2 - f2Dh >, (4.20) 


the inhnitesimal canonical transformation X^ generated by H is tangential both 
to Eh and Eh, and each integration curve of Xh on Eh H Eh yields a solution 
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to the Einstein equations. However, this correspondence is not one-to-one because 
the same spacetime allows an inhnite number of different slicings with the same N 
and u. Further for a different choice of the lapse function N and the shift vector u 
yields a different curve in E for the same spacetime. 

This ill correspondence between the spacetime solutions to the Einstein equa¬ 
tions and the curves in the phase space, which arises due to the general covariance 
of general relativity, can be made well-dehned by considering the subspace spanned 
by the integration curves in stead of each curve. To see this, let us denote the set 
of constraints symbolically by ha and the corresponding inhnitesimal canonical 
transformations by Ya = Xh^- Then from the first class nature of ha shown above, 
{Ya} yields an involutive system on fl T^h- 

(4-21) 

7 

where is a set of functions on F. Hence we obtain a foliation of the constraint 
submanifold Sj) fl = UaCa where each leaf Cx is a connected component of the 
integration submanifolds. 

For an arbitrary non-degenerate curve 7 contained in a leaf Cx, its tangent 
vector X is written in terms of some set of functions TV" as X = J2a N°‘Ya because 
{Ya} spans the tangent space of Cx at any point. Hence it is an integration curve 
of the hamiltonian flow for the hamiltonian H = J2aX°‘ha, and corresponds to 
some spacetime solution (M, g) to the Einstein equations. Further, if two curves 
7 i and 72 are contained in the same leaf and intersect with each other at a point u, 
they correspond to solutions to the Einstein equations with the same initial data 
for different lapse functions and shift vectors. Hence from the uniqueness of the 
initial value problem for the Einstein equations the corresponding spacetime solu¬ 
tions {M,gi) and {M,g 2 ) are isometric if they are maximally extended. The same 
conclusion holds even if these two curves do not intersect. For one can End an¬ 
other curve 73 in the same leaf which intersects both with 71 and 72 , which implies 
that {M,g 3 ) = {M,gi) and {M,g 3 ) = {M,g 2 ), hence {M,gi) = {M,g 2 ). Therefore 
all the curves contained in the same leaf corresponds to a unique 4-dimensional 
diffeomorphism class of the spacetime solutions to the Einstein equations. 

We can further show that this correspondence is one-to-one. Take two curves 
7 i C Cl and 72 C C2 and suppose that the corresponding spacetime solutions 
{M,gi) and {M,g 2 ) are isometric. Then there exist isometric diffeomorphisms to 
a spacetime {M,g), <hi : {M,gi) {M,g) and <h 2 : {M,g 2 ) {M,g). Let Si = 

X S) and S 2 = $2 ({^ 2 } x S) be two space-like constant-time hypersurfaces 
in M, and choose two families of time slicings in M such that they contain a 
common time slice, one of them contains Si and the other 5 * 2 . Further let the two 
curves in the phase space determined by these slicings be 73 and 74 . Then from the 
construction 73 ft 71 7 ^ 0, 74 fl 72 7 ^ 0 and 73 H 74 7 ^ 0. This implies that there is a 
curve which connects a point in Ci and a point 6 * 2 . Hence from the connectedness 
of each leaf Ci and C 2 must coincide with each other. 
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Thus we have found that the connected components of the integration man¬ 
ifolds of the involutive system {Ya} are in one-to-one correspondence with the 
4-dimensional diffeomorphism classes of the spacetime solutions to the Einstein 
equations. In other words the classical dynamics of general relativity is completely 
determined by the foliation of the constraint submanifold in terms of the inhnites- 
imal canonical transformations generated by the constraints. We can discard the 
lapse function and the shift vector, or the corresponding hamiltonian. For this 
reason we can simply say that the canonical theory of general relativity is given by 
a totally constrained system (T*A4{S),uj, {ho, hj.}). We will call each leaf of the 
foliation a causal submanifold. 


4.3 Elimination of the diffeomorphism constraints 


As the structure constants of the Poisson brackets among hn are genuinely con¬ 
stant from Eq.( [4.18| ), the corresponding inhnitesimal canonical transformations are 
involutive on the whole phase space and generate the action of the diffeomorphism 
group of S', Diffo(S') where the subscript 0 denotes the connected component con¬ 
taining the unit element. Since all the measurable quantities are invariant under 
these transformations, it is desirable to eliminate this kinematical gauge symmetry 
from the canonical theory, especially when one consider the quantization of the 
theory. Now we will show that the classical dynamics has the same structure as 
above even after the elimination of this gauge freedom. 

First of all note that for a jF(S')-valued functional iV on P which transforms 
covariantly under Diffo(S') as 


N{a^u) = a,{N{u)) Va e Diffo(^), Vu e P, (4.22) 

< hx, iV > is invariant under Diffo(S') as a functional on P from 

h_|_, N >, < hoi ^ ^ j" ^xh±i -Y > -|- < h_|_, £;xN >= 0. (4.23) 

Similarly for a r¥(S')-valued functional h on P which is Diffo(S')-covariant, < Iid, v > 
is invariant under Diffo(S'). Further by inspecting the argument on the correspon¬ 
dence between a curve in a causal submanifold and the hamiltonian flow generated 
by the hamiltonian H =< h±,N > + < hD,i^ > one easily sees that N and 
z/ can be replaced by some appropriate Diffo(S')-covariant functionals N and h. 
Hence the connected integration surfaces of the involutive system generated by the 
Diffo(S')-invariant functionals < h±,N > and < h > give the same foliation as 
that given by < h^, u > and < h±, N >. 

Further if < h£),f 0 at a point u G P for some f G A’(S'), there exists a 

functional : P —> A’(S') such that < liD^ i 0 the same point u. Thus Sj) 
can be redehned as 


Sd = {"W G P I < ho, p > iu) = OVz/ : P —r¥(S'); Diffo(S')—covariant} . (4.24) 
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Similarly Tih can be expressed as 

e r I <h^,N>{u)= OViV : r Diffo(^)-covariant} . (4.25) 

These arguments indicate that the original canonical system can be naturally 
projected on T /Diffo(S'). To conhrm this, let us denote all the functions on T which 
are invariant under Diffo(S') by 

^inv := {/ G I {/, < ho,^ >} = ove e ^(5)} . (4.26) 

Then it is easily shown that is closed with respect to the Poisson algebra and 

where 

l/ls, = 0}. (4.27) 

Further since {/, < hr,,^ >} = 0 implies that Xf is tangential to T^d, {f,g} = 
—Xfg vanishes on T^d for / G JFinv and g G Xd- Hence {XirwT^n} = ^d- This 
implies that Xd forms an ideal of jFjnv and the Poisson bracket in jFjnv naturally 
induces a Poisson bracket in A-mv '■= X-mvlX d- Each element of Aim is just a set 
of functions in Xim which coincide with each other on 

Let TT : P —P/Diffo(S') be the natural projection and put Pinv := Then 

from the arguments above Pinv is characterized as 

Tinv = {"W G P/Diffo(S') I < h^), v > (u) = 0Vz> ; P —> ?T(S'); Diffo(5')—covariant} , 

(4.28) 

and Aim is naturally identihed with jF(Pinv). Further the constraint ho is trivial¬ 
ized on Finv and the causal submanifolds in fl is bijectively mapped to the 
causal submanifolds in ti{T,h) H Finv determined from < h_L, iV > G Finv 
Let ha G Aim be a generating set of all the functions of the form < hj_, 
such that 

i) For any Diffo(S')-covariant functional N : T ^{S) there exists a set of 
elements A“ G ^inv such that < h_i, N >= J2a 

ii) 7r(S/^) n Finv = {m G F^v | ha{u) = OVa}, 
hi) {ha.hn} 

Further let us denote fl Finv by the same symbol T^h- Then the arguments 

so far shows that the canonical dynamics of general relativity is described by the 
totally constrained system (Finv, cuinv, {ha}) and the causal submanifolds in T^h is 
one-to-one correspondence with the Diffo(M)-class of the spacetime solutions to 
the Einstein equations. 
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Figure 1: 


4.4 Cotangent-bundle structure of Finv 

In the last statement of the previous subsection (Ujnv is understood to be the sym- 
plectic form corresponding to the Poisson brackets in Aim- Hence in order to make 
the statement rigorous it should be shown that Pinv has a manifold structure and 
the required symplectic form exits. Now we prove these facts by showing that Pinv 
can be identihed with T*(A^(S')/Diffo(S')) and cUmv coincides with the canonical 
symplectic form corresponding to the cotangent bundle. 

Let 01 : Ed —> M(S) be the restriction of the natural projection from T*(Af (S')) 
to Af(S). Then, since the diffeomorphism constraint implies that p vanishes on 
the subspace of T*(A^(S)) spanned by the tangent vectors to the Diffo(S)-orbits 
from Eq.( |4.7| ), 0i is surjective and induces a surjective mapping 02 : Tiny —> 
Af (S)/Diffo(S) such that 7r0i = 027r(See the diagram in FigJ^). Let G Tg(A4{S)) 
be a vector tangent to a Diffo(S)-orbit passing through q G M(S). Then, since it 
is written in terms of a vector held X G A’(S) as pg = it follows from Eq. ([4.7|) 
that for (p, g) G 


(P, ?) ijlq) =< P, Pq >=< P, ^xq >= 0. (4.29) 

Further for a G Diffo(S) and v G Tq{M.{S)) 

a*(p,g)(a*u) =< a^p,a^v >=<p,v >= (p,g)(u), (4.30) 

from the diffeomorphism invariance of < p,v >. Hence there is an injection j : 
Finv —^ ^*(A^(S)/Diffo(S)) such that for v G Tq{A4{S)) 

j7r(p,g)(7ru) =< p,u > . (4.31) 

It is easily checked that for the natural projection 03 : T*(Al(S')/Diffo(S')) — >• 
A4{S) /Diffo(S'), 02 = (psj holds and that j is surjective. Thus Finv can be identihed 
with T*(Al(^)/Diho(^)). 

Next let us show that the symplectic form cuinv induced from the cotangent 
bundle structure of F^v is equivalent to the Poisson brackets in Amv — A(Pinv) 
derived from the symplectic form uj in P. From now on we identify Finv with 
T*(AI(S')/Diho(S')) and write jir simply as tt. 
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We first show that each element of ^inv uniquely determines a vector held on 
T*(Wl(S')/Diffo(S')), which will turn out to be an inhnitesimal canonical transfor¬ 
mation on Finv Let f¥inv be 


A’inv = {X e X{Ed) \a,X = XWae Diffo(^)} , (4.32) 

and for [/] G Anv, let r¥[j] be 

= {Xe A’inv I u;{x, Z) = -Zf VZ e . ( 4 . 33 ) 

Then for each X G obviously dehnes a unique vector held on T*(Af (S')/Diho(S')). 

Further let 6* and 6*inv be the canonical 1-forms on T*{Ai{S)) and T*(Af (S')/Diho(S')), 
respectively. Then from the commutativity of the diagram in Fig.[^ and Eq. ([4.31| ) 
it follows that (6'inv)7r(p,g)(7r*X) = 9(^p^q){X) for any X G X-mv, i-e., 

e{X) = 7rW,^,{X) MX G ^inv (4.34) 

Hence for AX = Xi — X 2 (Xi,X 2 G ^[f]) and Z G ^mv, from 

0 = cu(Ax, z) = de{Ax, z) = Ax{e{z)) - z{e{Ax)) - 0([ax, z]), ( 4 . 35 ) 

it follows that 

0 = 7r,(AX)(0i„.(7r,Z)) - 7r,Z(0i„.(7r,X)) - 0i„v(7r4AX, Z]) 

= d6*inv(vr*AX, 7r*Z) = a;inv(7r*AX, 7r*Z). (4.36) 


As TT^Z can be any vector held on T*(Af (S')/Diho(S')), this equation implies that 
vr^AX = 0, i.e., tt^Xi = 7 r*X2. Hence [/] G ^inv determines a unique vector held 
on T*(Xl(S')/Diho(S')). We denote this vector held by X[j]. 

Next we show that for X G r¥[j] and Y G r¥[g] a;(X, Y) = —{/, g}. For Xi, X 2 G 
r¥[/] and yi,T 2 £ ^[g] from the dehnition (^4.33|) it follows that 

o^(X2, Y 2 ) - a;(Xi, W) = a;(X2 - Xi, Y^) + o^(Xi, ^2 - W) = 0. (4.37) 


Hence a;(X, F) depends only on [/] and [g]. Obviously Xj G X[f] and Xg G ?i[g\. 
Therefore 

{/, g} = -uiXf, Xg) = -uiX, Y). (4.38) 

With the help of the equations derived so far for X G r¥[j-] and Y & ^[g\ we 
obtain 


[{f,g}] =-[u;{X,Y)] = -[d9iX,Y)] 

= [Yi9 iX))-Yi 9 iX))- 9i[X,Y])] 

= 7r*F(0inv(7r*X)) - 7r,F(0i„,(7r,X)) - 9inA^,[X, F]) 
d0[yyy(^7r^X^ ) 

= {|/].l9]}in. (4.39) 
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This shows that the Poisson brackets induced from oj coincides with that defined 
by a;inv 

Note here that the arguments so far are not mathematically rigorous because 
A1(S')/Diffo(S') has conical singularities at metrics with Killing vectors@. Though 
these singularities may have physical importance in quantization, we will not go 
into this problem in this paper. 

We can go further and eliminate all the hamiltonian constraints to get the fully 
reduced phase space with a symplectic structure which represents the true physical 
degrees of freedom as done by Fischer and Marsden. However, we shall not follow 
this line because we will then lose the dynamics. 

5 General Totally Constrained Systems 

Now we discuss the dynamics of a generic totally constrained system. Here a totally 
constrained system is defined as a triplet of a phase space, a symplectic form and a 
set of constraint functions, (r,a;, {h„}). For a technical reason we assume that the 
phase space is 2n-dimensional smooth manifold with finite n. Further we assume 
that the constraints are of first class with the Poisson brackets given by 

{ha,h0} = (5.1) 

7 

where are functions on F. 

On the basis of the arguments in the previous section we understand that the 
physical evolution of the system is one-to-one correspondence with each leaf of the 
foliation determined by the involutive system of the infinitesimal canonical trans¬ 
formations Ya = on the constraint submanifold T,h = {m G F | ha{u) = 0 Va}. 
We call each leaf a causal submanifold as so far. As is clear from the arguments 
in the previous section, this interpretation is equivalent to regard that two solu¬ 
tions to the canonical equation of motion for the hamiltonian H = X°ha with 
arbitrary functions A" represent the same physical evolution if they intersect with 
each other in F. 

This is a natural generalization of the argument on the dynamics of a single 
totally constrained system with one constraint in §3. Now we extend this general¬ 
ization to the statistical dynamics of an ensemble. 

5.1 Relative distribution function 

From this interpretation of dynamics of a single system and the argument in §3 
it is natural to introduce the relative distribution function p on F to describe an 
ensemble, which vanishes outside the constraint submanifold and is constant on 
each causal submanifold: 


YaP = 0 Va, 


( 6 . 2 ) 
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haP = 0 Va. 


( 5 . 3 ) 


Let us define an acausal submanifold as a submanifold of L which intersects with 
causal submanifolds transversally. Then for any acausal submanifold S and for 
any distribution p|s on S a solution to these equations which coincides with p|s 
on S is unique, if it exists, on the causal development of S dehned by 

£>(E):= U C (5.4) 

CnS7^0 

where C runs over causal submanifolds. However, such solution may not exists in 
general. In fact the following theorem holds. 

Theorem 5.1 In order that there exists a solution to Eqs. for arbitrary 

initial data on any acausal submanifold, the following condition should be satisfied: 

- 0 - 

7 

This condition is satisfied if and only if there exists a function f ^ 0 such that for 
c '^/3 corresponding to the constraints h'^ = fh^ 

c'a := Eca/3 ~ 0 
P 


holds. 


Proof 


Since p is a distribution, to be exact, Eqs. (|5.2|) - (|5.3| ) are expressed as 


< haP, 0 >:=< p, ha4> >= 0, 

< YaP, (j) >:=< p, -Yacf >= 0, 


where (f is an arbitrary smooth function with a compact support on T. However, 
since the commutators among Y^’s are given by 

7 

the consistency condition yields 

0 =< p, |y„, Yp],!, >=< p, Y.(y;clp)4> > . 

7 

Hence, noting the relation Y^c^^p = —{h^, 'w® obtain the first condition in the 

theorem. 
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In order to show the latter half of the theorem, hrst note that the Jacobi identity 
for the Poisson brackets among ha yields 

YaCp - YfSCa + Y.CI0- , 

7 7 

where Cq, := From this it immediately follows that the hrst condition of 

the theorem holds if ~ 0. 

On the other hand for h'^ = fha Ca changes as 

Ca = Ca - (m- l)Ya\nf, 

where m is the number of the constraints. Hence the second condition of the 
theorem is satished if / is a solution to the equation 

Yar-^ ~ 


However, if the first condition of the theorem is satisfied, we obtain 




-Ec. 


7 

a0' 


This is nothing but the consistency condition for the hrst-order differential equation 
system for above. Hence the hrst condition of the theorem implies the second. 


Note that for a matrix function A = (A^) on F with det A 7^ 0, the totally 
constrained system with the constraints h'a = J 2 p-^a ^0 is equivalent to the original 
system. Hence the precise meaning of the requirement of the theorem is that 


Cap} can be put to zero by such a transformation and that Eqs. (|5.2|) -( 
are consistent only for such choice of the constraints. This result is interesting in 
relation to the quantization of the totally constrained system because this condition 
implies that the operators corresponding to and ha should commute in a weak 
sense. 

On the basis of this theorem we assume that = 0 from now on. Under this 
condition if we put 

P = (5.5) 

a 

Eq. (^.2|) is automatically satished and Eq. (|5.3|) reduces to the equation for the 
function pQ, 

YaPo ~ 0, (5.6) 

because 


= U n + T. <r,h-,S'{h) n 

0 0n 

= (XaPo - CaPo)X{h{h^). 


(5.7) 
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5.2 Statistical dynamics in terms of conservative measure 
on acausal submanifolds 

Now we show that we can formulate the statistical dynamics for an ensemble of 
the totally constrained system without reducing it to a system without constraint. 
The basic idea is the same as that used in §2 and §3 for totally constrained systems 
with a single constraint. 

First note that for any pair of acausal submanifolds Si and S 2 the foliation 
of the constraint submanifold by the constraints uniquely determines the causal 
mapping 

9 Si n d(Yj 2 ') n Yjjj —> S 2 n z^(Si) n s^^^, (5-8) 

where T*(Si) and T*(S 2 ) are causal developments of Si and S 2 , respectively. We 
extend this causal mapping to a neighborhood of by considering a foliation of 
the tubular neighborhood such that the intersection of each leaf with Sj^ coincides 
with the foliation of Sj^ by Tq- If a measure pso on an acausal submanifold Sq with 
its support contained in Sq fl T,h is given, this causal mapping uniquely determines 
a measure /is on S fl T*(So) with its support contained in such that for any 
constant of motion, i.e., a function / G JF(r) which is constant along each leaf, 

/ fdfiT. = [ /d/iSo, (5-9) 

JS JSo 

if supp / n So C dom 9. As in the single constraint systems, this measure can be 
expressed locally in terms of p, and fl. 

Theorem 5.2 If Ca = 0, for any distribution p satisfying Eqs. 
sure 

d/is = p\Iyi ■ ■ -dy^flls 
is conserved by causal mappings where m is the number of independent constraints. 

Proof 

Let us denote ■ ■ ■ /y^ (pfl) simply as y. Then from the identities 

d o Jx + /x o d = &X-, 

[^x, dy] = 

we obtain 

m 

dx = -lY.dly, ■ ■ ■ Irjpn) + Y, cl„(-i)“ n lYgipn). 

a=l (3^0. 

In a similar way we obtain 

(_l)7-l/y^ . . . . . . /y^(pfl) = (-1)7 ■ ■ ■ (pfl) 

m 

+(-ir' E . 

q := 7+1 /3^'y 


(5.5) the mea- 
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Since d{pQ) = 0 and £^y^(pf2) = 0, summation of these equations over 7 = 1 ,..., m 
yields 

m 

<ix = E(-ir'''=7EAi(pn) = o. 

7=1 / 3^7 

Hence for a pair of acausal submanifolds Si and S 2 such that -D(Si) = D{T, 2 ), 
from Stokes’ theorem on {2n — m + l)-dimensional submanifold N such that dN = 
Si U S 2 U S' and S' is parallel to the leaves, we obtain 

/ \x\- [ \x\ = ±f XT f X = ± f dx = 0. 

Jt.2 -'Si JdN Jt,' Jn 


In realistic situations each acausal submanifold is specihed by a set of indepen¬ 
dent m functions 0a such that Ya(j)p = {0/?, ha] is a regular matrix as, say. 


h — ■ ■ ■ — 0m — 0. 


(5.10) 


Let us call such a set of functions instant functions. Then the measure given in 
the previous theorem is expressed in terms of these instant functions as follows: 

Theorem 5.3 // 0i, • • •, 0m are instant functions for an acausal submanifold S, 
for any / G JF(r) the following equality holds: 

[ f p\Iyi ■ ■ ■ = f /p|det{ha,0/3}in'^(</>7)l^l- 

Proof 

From 

d(f)i A ■ ■ ■ A dcfrn A Iyi ■ ■ ■ 

A ■ ■ ■ A df)a^ My,--- Iy^^^ 


ml 


ai---an 

m 


■ 1 ) 


m 


m\ 

(- 1 ) 

{m — 1)! 
(- 1 )”^ 


H e“iA ■ ■ ■ A ci0a^ A 1^2 ■ ■ ■ Iy^^) 


e“l'"“’"(/yjd0ajci0a2 A ■ ■ ■ A d0a^ A Iy^ ' ' ' 


(m- 1)! 


Y e“i'"“’"{/ii,0ajd0a2 A • ■ ■ A d4>ara A Iy2 ' ■ ■ 


we obtain 


Qf ■■■CXrn 


I[d{(j)a)Y[d(l)p\lY, ■ ■ ■ lYm^ = I det{ha,0p}| 



By multiplying fp on the both sides of this equation and integrating over F, we 
obtain the equation in the theorem. I 

From these theorems we can formulate the statistical dynamics of an ensemble of 
the totally constrained system with multiple constraints in the following way. First, 
from the data set obtained by measurements, pick up a set of instant functions 0(i)a 
which take a common set of values in the data set. For simplicity assume that these 
values are all zero, and let the corresponding acausal submanifold in F be Si, and 
dehne the measure du by 

diy:= c\Iyi ■ ■ ■ lYm^l (5-11) 

where c is some hxed positive constant of motion. Then the other data uniquely 
determines the distribution p on Si through the formula 


< / >si= / fpdy. 

JEi 


(5.12) 


Extend this distribution p over D(T,i) by the evolution equations Eqs.( b-2|) -( b.3|) . 
Then for another set of instant functions (p( 2 )a corresponding to an acausal sub¬ 
manifold S 2 , the expectation value of a function / G JF(F) on that submanifold is 
given by 


< / >E2= / fpdu, 


IT,2 


(5.13) 


if supp / G /l(Si). Of course we do not need the explicit knowledge on the acausal 
submanifolds, because from Theorem ^.3| the expectation values are written as 
an integration over F in terms of measures expressed by the constraints and the 
instant functions. 

Like the case of a single constraint system we can dehne that a set of functions 
4>a are good time variables if the natural measure — '7‘a)|0| on a set of 

acausal submanifolds 

S(r„) := {u G F I = r„} (5.14) 


coincides with the conserved measure. This condition is expressed as 


c| det{h„,<l>p}\ K -Am), 


(5.16) 


where c is some positive constant of motion and / is some function of m variables. 


6 Discussion 

In this paper we have shown that the dynamics of a classical totally constrained 
system can be consistently formulated without reducing it to an unconstrained 
system by solving the constraints or referring to a special time variable. The basis 
idea has been to consider the relative distribution function which is constant on 
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each leaf of the foliation defined by the inhnitesimal canonical transformations gen¬ 
erated by the constraint functions, and to normalize it on an acausal submanifold 
which is transversal to the foliation in terms of the conservative measure. 

The fact that we can formulate the classical statistical dynamics of a totally 
constrained system without referring a special time variable is very important for 
considering a quantum theory of the totally constrained system because a quantum 
theory has a similar structure to the classical statistical dynamics in general. In fact 
in the next paper we will show that by introducing a similar foliation structure into 
a state space of quantum theory and by considering a relative probability amplitude 
in stead of the relative distribution we can construct a consistent formulation of the 
quantum dynamics of a totally constrained system without referring to a special 
time variable under some restrictions. 

Though the main purpose of the present paper has been to give a basic motiva¬ 
tion for the quantum formulation developed in the subsequent papers, the results 
obtained in the paper may be interesting by themselves. In particular the fact that 
the conservative measure can be written only by the canonical volume form and 
the constraint functions even for multi-constrained systems seems to be useful in 
the arguments of the probability distribution of the initial condition of the universe 
in the classical framework and stochastic treatment of general relativity. 

Of course the expression for the measure given in this paper cannot be applied 
to general relativity directly because we have only considered systems with hnite 
degrees of freedom. However, it seems possible to extend the formulation to gen¬ 
eral relativity by taking an appropriate limit. To examine this limiting procedure 
explicitly in some simple situations such as the perturbation theory of general rela¬ 
tivity on cosmological background spacetimes and spherical black hole spacetimes 
with scalar helds will be interesting. 
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